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1 Introduction 



With the evidence in favor of neutrino masses and mixing at a convincing level now, attempts 
to seek the role of physics beyond the standard model in the neutrino sector are acquiring 
enhanced degrees of urgency. As it is, the lack of naturalness of the mass of the Higgs 
boson in the standard model is a strong pointer towards new physics around the TeV scale. 
Since tiny masses can be explained by appealing to energy scales much higher than the 
electroweak scale (for example in the seesaw mechanism), it is appropriate to link neutrino 
mass generation to new physics options at the TeV scale or above. 

Supersymmetry (SUSY) is a frequently explored possibility for new TeV-scale physics. 
Its capability for solving the naturalness problem being an accepted fact, serious efforts are 
on at accelerators to see signals of SUSY, broken with an intra-supermultiplet mass splitting 
O(TeV). Does SUSY play a role in providing the requisite new physics component in the 
masses and mixing pattern of neutrinos? This is the question that we would like to address. 
Often one has to go beyond the minimal SUSY standard model (MSSM) in order to find 
satisfactory mechanisms which can achieve this end. Though a fair amount of work has been 
done in this area pQ , one is yet to have a satisfactory answer to the following central question 
related to the neutrino sector. Why is the mixing pattern of neutrinos, with two large and 
one small mixing angles, so drastically different from that in the quark sector, where the 
mixing between the three generations can be cryptically described by the progression - small, 
smaller, smallest? 

In this paper we approach the above problem with the idea that the difference between 
the two mixing patterns arises from some aspects of the SUSY model which are specific to 
neutrinos and with no counterparts in the quark sector. For this, we make use of nonrenor- 
malizable terms arising from high-scale physics. Such terms, coupling some hidden sector 
(gauge singlet) chiral superfields to the MSSM ones, are suppressed by the Planck mass 
Mp or some power of it. If these terms violate lepton number , they can lead to Majorana 
masses for neutrinos. When, in addition, there are superfields containing right-chiral neu- 
trinos, contributions to the neutrino mass matrix can come not only from the well-known 
seesaw mechanism but also radiatively via one-loop diagrams containing right chiral sneu- 
trinos. Though both these contributions have been included in earlier works [2J El El a 
clear explanation of the different character of mixing for neutrinos vis-a-vis quarks has been 
lacking without the imposition of some additional restriction on the low-energy theory. The 
Froggatt-Nielsen mechanism [3] is an example of such additional theoretical inputs. The 
literature, of course, is rich with uses of various other symmetries [6 , as well as of 'anarchy' 
in the neutrino mass matrix jjj. Drawing inspiration from all these approaches, we suggest 
an alternative justification of bilarge neutrino mixing by postulating an array of gauge sin- 
glet chiral superfields with flavor- dependent nonrenormalizable couplings to neutrinos. The 
specific superpotential that yields the desired results is formulated, and consistency with 
the observed suppression of flavour-changing neutral currents (FCNC) in the lepton sector is 
used as a constraint. One further has radiative contributions, pertaining only to the neutrino 
mass matrix. These are due to the fact that the right chiral sneutrinos may acquire gauge 
singlet AL = 2 mass terms. The scenario for neutrinos then immediately becomes quite 
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distinct from that in the quark sector. 

Using the standard seesaw as well as the above-mentioned radiative contributions, we have 
examined whether high-scale parameters, such as the vacuum expectation values (VEV) of 
scalar as well as of auxiliary components of the gauge singlet chiral superfields, crucial to 
this mechanism, are in otherwise acceptable ranges of values. For instance, the Higgsino 
mass (//) parameter needs to be around the weak scale for the desired implementation of 
the spontaneous breakdown of electroweak symmetry. Such analyses are carried out for 
the three alternative possibilitiess of the neutrino mass spectrum allowed by the neutrino 
oscillation data, namely, normal hierarchy, inverted hierarchy and degenerate neutrinos. The 
simultaneous importance of the radiative as well as the seesaw contributions enables us to 
acquire in all the three scenarios substantial regions (of different extent in each case) in the 
parameter space of our model that correspond to acceptable solutions. 

In section 2 we describe the (by now well-known) structure of the neutrino mass matrix 
M v for bilarge mixing. The SUSY model is constructed and the elements of M v are consis- 
tently generated in section 3; we also show at the end of this section how FCNC processes, 
induced at one loop, are suppressed. The SUSY parameter space, answering to each of the 
specific scenarios of normal hierarchy, inverted hierarchy and degenerate neutrinos, is ana- 
lyzed in section 4. We comment on some related possibilities in section 5. Section 6 contains 
our summary and conclusions. 



2 Facts about neutrinos 



There is experimental evidence now that neutrinos have tiny masses. We shall work within 
the scheme of three light active neutrinos, not including the possibility of an additional light 
sterile one suggested by the LSND data till results from the ongoing mini-Boone experiment 
settle the issue. We shall further assume CPT conservation. While there is an upper bound 
|Hj of ~ 1 eV on the sum of neutrino mass eigenvlaues from cosmology, the lack of observation 
of neutrinoless double beta decay implies an upper bound of ~ 0.3 eV on the absolute value 
of the the ii-element of the neutrino mass matrix |9 3 . On the other hand, the accumulating 
data from solar, atmospheric, accelerator and reactor neutrino experiments persistently point 
[TO] towards neutrino oscillations. These data identify favored regions of small but distinct 
mass-squared separation of the three different physical neutrino states. At the same time, 
they also indicate that the mixing between the second and the third families is near maximal, 
that between the first and the second is large, while the one between the first and the third 
families is restricted to a small angle. In perfect analogy with the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix in the quark sector, the three-flavour neutrino mixing matrix can 
be parameterized as 
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In Eq. (1) em- 
phases have been neglected here). We work in the basis where the charged lepton mass matrix 



cos 8ij, Sij = sin 6ij, i, j being family indices which run from 1 to 3 (Majorana 
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is diagonalized. While solar (reactor) neutrino (antineutrino) studies suggest that 6*12 — 32° 
[TT| IT2"] , the atmospheric neutrino deficit needs #23 to be ~ 45° ^3] and data from reactors 
require that #13 < 13° pi] . Thus a pattern of bilarge mixing emerges. 

The above pattern allows one to construct a candidate neutrino mass matrix in terms of 
the mass eigenvalues mi,m 2 ,m 3 . To start with, let us take #23 = f and 6*13 = approxi- 
mately, keeping # 12 free to be large. Then the corresponding transformation matrix can be 
written as 



\ 



e 
V2 



1 

73 J 



(2) 



V2 

s _ c 

where s = sin# 12 and c = cos6 l 1 2. The neutrino Majorana mass matrix in the flavor basis 
can now be obtained by transforming the diagonal matrix with the above: 

\ 
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m 2 c 2 + m 3 ) \(mis 2 + m 2 c 2 + m 3 ) y 



Thus we see that the requirement of bilarge mixing commits one to a particular structure 
of the mass matrix where, of course, the relative magnitudes of the entries depend on the 
eigenvalues. For more precise information one has to take up the specific scenario of nor- 
mal/inverted hierarchy or that of degenerate neutrinos. In our study, we attempt to link the 
diagonal and off-diagonal mass terms of M v to the parameters of the SUSY model at high 
scale and see what the different scenarios tell us about the model parameters themselves. 



3 The SUSY model and neutrino masses 
3.1 Required features of the model 

The model that we adopt is motivated by a number of recent works El H|- I n |H| for 
example, a minimal extension of the MSSM, including a right-handed neutrino, is used. 
There the terms of the effective Lagrangian responsible for neutrino masses are 

Ceff = ^- ([X*NN] D + [XLNH U ] F ) + h.c, (4) 

where Mp is Planck scale and coupling coefficients of order unity have been suppressed. In 
Eq. (4), the chiral field X can acquire both SUSY violating and SUSY conserving VEVs [To'] . 
The above terms can be responsible for seesaw masses of order for the neutrinos (recalling 
the need to have Fx ~ (x) 2 ~ Mp(H u ), to ensure the generation of other superparticle 
masses in the TeV range) . In addition, there can be radiative contributions to the mass 
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matrix from AL = 2 sneutrino mass terms after SUSY breaking. These can arise with 
the help of a term j^-[X^XX'NN]d, making contributions that can dominate over seesaw 
masses in certain regions of the parameter space. 

We aim to explain the bilarge mixing pattern of neutrinos by extending such a model. As 
mentioned earlier, the basic philosophy is to envision some feature of neutrinos, which has no 
counterpart in the quark sector, as being responsible for the observed bilarge mixing. There 
are two features of this kind in such a model: (a) the right-chiral neutrino sector and (b) the 
corresponding right-chiral sneutrino sector, with provisions of AL = 2 terms in each. In our 
approach, each of these sectors is attributed with a 3 x 3 mass matrix structure which plays 
a crucial role in the contribution to the radiative as well as seesaw mass terms. Moreover, 
we postulate an array of gauge singlet chiral superfields Xij. Following these propositions, 
Eq.(4) is generalized to 



where there is a summation over flavor (i.e. generation) indices i, j. The choice of the above 
Lagrangian can be motivated by a global symmetry GpxG j2]. The factor G helps in solving 
the /i-problem, and keeps the spontaneous SUSY breaking scale \JF X low enough for the 
superparticle spectrum in the observable sector to be around TeV energies. The summation 
over family indices in Eq. (5) can be justified by the global symmetry Gp- This essentially 
means that the hidden sector chiral superfields X^ interact with those of the visible sector 
with such a global symmetry and that the low-energy flavour structure is the artifact of 
such interactions. In case the X^s acquire SUSY violating VEVs, then different soft SUSY 
breaking terms will arise from the nonrenormalizable interactions shown in Eq. (5) and 
other higher dimension terms compatible with all symmetries of the theory. These arise in 
addition to the soft terms that have analogues in the squark sector. 

Schematic expressions can be written for the neutrino and sneutrino mass terms, thus 
obtained, and for the soft SUSY breaking A-terms as well as for the corresponding terms in 
the SUSY Lagrangian leading to them. They are obtainable from the following realizations. 
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mj(Fx kl )(x* l3 ) + (F^)(F Xkl )(x* k ) + (F Xki )(F x J(x* l3 ) + 
(F Xl .){F Xik )(x* kl ) + (F Xik )(F xlj ){xti) + (F Xkl )(F Xlj )(x* k )}, 



6Mj> 
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where Mew — v/\^2 = 17 '4 Gev. The above expressions are up to unknown multiplicative 
factors occurring in the SUSY Lagrangian. It is also assumed that the Dirac mass matrix, 
generated by canonical Yukawa couplings (as in the quark sector) arising from renormalizable 
terms in the superpotential, has very small off-diagonal elements. Also, in addition to the 
AL = 2 mass terms shown above, there may be L-conserving mass terms for right-chiral 
sneutrinos result of soft SUSY breaking. 

Nondiagonal A-terms can potentially contribute to FCNC processes such as fi — > and 
hence need to be suppressed. Therefore, we wish to have a structure where Fx i vanishes for 
i ^ j. On the other hand, contributions to off-diagonal terms in the neutrino mass matrix 
are essential for bilarge mixing. Such terms will be made to arise from seesaw as well as 
radiative processes. As we shall show below, both of these are driven by the VEVs of the 
scalar components of the X-superfields. We must therefore have nonzero (xij) for i ^ j. 

Thus we require nonrenormalizable terms in the superpotential involving the chiral super- 
fields X^, with a rather interesting complementarity between the diagonal and nondiagonal 
elements of the array. The diagonal ones can have nonvanishing F-term VEVs and thus can 
generate a diagonal A-matrix, whereas the off-diagonal elements must have vanishing F-term 
VEVs, though the corresponding scalar VEVs must be nonvanishing. A superpotential, in 
which the above characteristics can be achieved, is presented below. 

3.2 The superpotential 

It was noticed in the previous subsection that bosonic components of the array of chiral 
superfields should acquire SUSY violating and SUSY conserving VEVs in a complemen- 
tary manner to be able to generate the required neutrino masses and mixing pattern. In 
order to achieve this, we first demonstrate a simple situation in which the auxiliary and the 
scalar components of a single chiral superfield acquire SUSY violating and SUSY conserving 
VEVs respectively. Thereafter we generalize this to an array of such superfields X^. 

Consider a set of hidden sector fields, for which the superpotential is of the form E] 

W = S(YY - n\) + Y 2 X' + Y 2 X 

+S'(XX - fj%) + X 2 Z + X 2 Z. (10) 

Here the i?-charge for the chiral field X is |, while the .R-charges of the remaining chiral 
fields can be chosen so that W has -R-charge 2. The explicit assignment of -R-charges will 
be shown after generalizing this to an array X^. The superpotential W, shown above, leads 
to a scalar potential V which has local minima. The position of the abosolute minimum 
depends on the parameters /ii and /i2. 

Case(l): If < | 1 , the true minimum occurs at (y) = (y) = (s) = (s') = (x f ) = 
(z) — (z) = and = = In this case we have (x) ^ but F x = 0. Thus 
there is a nonzero scalar VEV, but SUSY breaking does not show up in the obserbable sector. 

Case(2): If |/i 2 | < |/xi|, the true minimum occurs at \(y) \ = \(y) \ = ^ and (x) = (x) = 
(s) = (s f ) = (x') = (z) = (z) = 0. In this case (x) = but F x ^ and SUSY is sponta- 
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Hidden sector 


Field 


V ; A,; X' ) ; ; )',: S; , S" • Z, ; Z; ; 
^Mj J J ij ^tj *->ij *->ij 

1 111 5 5 o o 8 4 

3 33 66 33 


R- charge 


Visible sector 


Field 


Qi L i Jji D i E i N i Ru Hd 

111111-1-1 
2 2 2226 


R- charge 



Table 1: R-charges of hidden and visible sector superfields. 

neously broken. 

Utilizing these lessons, we make a straightforward genralization of the above superpoten- 
tial to include an array of chiral superfields X^: 

w=jz w tJ , (ii) 

where 

= S ij (Y ij Y ij - //?•) + Y^X'^ + Y^Xij + 

S'^jXtj - ,,';) + X%Z {j + X%Z ir (12) 

The i?-charges of various superfields in this scheme are shown in Table 1. 

With the above R-charge assignments, the nonrenormalizable interactions relevant to us 
are all R- invariant. Moreover, the /^-parameter is obtained in the desired range from the 
term jL £ j(S% + S^)H,H d d*e, yielding 

(13) 

where F$ includes a sum over indices. Thus this scenario has the additional virtue of ex- 
plaining the value of fi around the electroweak scale, so long as the hidden sector F-terms 
can be justified to be at an intermediate scale ~ 10 11 GeV. 

The minima of the scalar potential arising from the above superpotential occur at (s^) = 
i s ij) = 0' = 0) ( z ij) = (zij) = an d further depend on the parameters /z^,/^-. We 
choose these parameters in the following way: 

1. For i — j, choose |/^-| < so that F Xij ^ and (x^) = 0. 

2. For % ^ j, choose \fj,ij\ < |//^-| so that F Xij = and (x^) ^ 0. 

The generation of off-diagonal entries in the neutrino Majorana mass matrix via nondiago- 
nal (xij) is ensured by this potential. On the other hand, we have secured a diagonal form 
for F Xij , thus suppressing contributions to FCNC processes from A-terms. This interest- 
ing complementarity is achieved rather naturally by postulating in the superpotential the 
presence of some mass parameters ^ and ■ and their relative hierarchies. Though these 
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parameters may all be broadly of the same order, the relative magnitudes of the primed and 
unprimed ones can naturally be quite different for different members of the array. The actual 
suppression of FCNC processes will be demonstrated in further detail in a later subsection. 



3.3 Neutrino mass matrix 

Schematic expressions for neutrino and sneutrino mass terms, induced in this scenario, have 
already been shown. Now we obtain the exact entries in the neutrino mass matrix. These 
will enable us to establish links between observable quantities and the parameters of the 
SUSY model. The superpotential yields Fx ijt = Fi^ij and (xij) = for i = j. For simplicity, 
we shall further assume that all VEVs are real and Fi = AMp for all i, thus reducing A to 
a single number. After SUSY and electroweak symmetry breaking, Eq. (5) then reduces to 

C eff = -AN l N* + A—uV - ^^N j + h.c, (14) 

where N stands for right-chiral neutrino fields (and not the corresponding superfields) . From 
Eq. (14) Dirac neutrino mass elements, as indicated already, are given by 

= j£T2 (15) 

while right-handed neutrino mass elements are given by 

Wij = 2AV (16) 

We can immediately deduce the seesaw masses from the above via the relation 

m s v = —mDTrip^m^). (17) 

In C e ff we could also use the term ^r[(Xj k XkiX^ + X ifc X^X/ J - + X^X^Xi^N 1 N j ] D as 
well as its hermitian conjugate which are consistent with all conserved quantum numbers. 
After SUSY breaking this term yields 

A 2 

£>hh = j^(xij)n l h J . (18) 
Consequently, the L-violating mass-squared terms for right-chiral sneutrinos UE>, become 

A?, = -A 2 ^. (19) 



The insertion of L-violating sneutrino masses allows the entry of radiative mass terms via 
the loop diagram shown in Fig. 1. The expression for the loop-induced contribution is 
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Figure 1: Radiative diagram at one loop level which generates Majorana neutrino mass. 



where fh is the SUSY breaking scale in the observable sector and is of the same order 
as the physical neutralino and sneutrino masses. The origin of the parameter A in this 
scenario has been shown in equation 8. It is sufficient for us to assume (x^'} — \P^j%l ^nd 
A = Ai = Fi/Mp = F/Mp for all i,j, which makes each of the above matrices symmetric 
in With this choice and after making the transformation v — > iv so as to change the 
overall sign of the neutrino mass term, the seesaw and radiative mass matrices respectively 
become 



m, 



AMpF 



( (xu) 2 + (x n ) 2 (xi 3 )(x 23 ) (x 12 )(x2 3 ) i 

(Xl 3 ){X23) (Xl 2 ) 2 + (X 2 3) 2 (Xl 2 )(Xl3) 

(xi 2 )(x 23 ) (xi 2 )(x 13 ) (x 13 ) 2 + (x 23 ) 2 J 



(21) 



384vr 2 Ml 2m 5 



( (x 12 ) (x V3 ) \ 

(x 12 ) (x 23 ) 
\ (xia) (x 23 ) J 



(22) 



In the above expressions we have used A = The uncertainty in caused by running 
clown to the electroweak scale can be absorbed in (xij) and F, since we are concerned 
with only the orders of magnitude of the latter. Also, the masses of all superparticles 
such as neutralinos and sneutrinos have been clubbed together as m here. With such an 
approximation already made, the effect of renormalization group evolution is not expected 
to make much difference. Finally, with the effects of both nonrenormalizable interactions 
and lepton- number violation, taken into account, our most general neutrino mass matrix is 



m v = ml + m r v . (23) 

Comparing the above mass matrix with Eq. (3), we are led to the following equations in 
the notation of section 2. 



m lC 2 + m 2 s 2 = ^-I((x 12 ) 2 + (x 13 ) 2 ) (24) 
1 v 2 1 

_( miS 2 + m2C 2 +m3) = ___(( Xl2 }2 + ^ 23 ^ (25) 
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= 4jrM Xl3}2 + {X23)2) (26) 

cs g 2 F 4 v 2 

-j=(-m 1 + m 2 ) = - OQ/ ,_ 9 ^75 o„--s ( a; i2) + TTT^^fa) (27) 

-^(m!-m 2 ) = - op " _ g ^ 5 l B <gi 3 ) + 777" ^faa) < g 23) (2? 

i(-m lS 2 -m 2 c 2 + m 3 ) = -7^775^(^3) + 777-75^2) (x 13 > (29) 

One set of consistent solutions to equations (25), (26) and (27), (28) is (xu) = —(^13). This 
reduces the above six equations to four which can be expressed as follows: 
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v 2 |(x 12 )| 2 t; 2 |(x 12 )| 2 v 2 |(x 23 )| 2 

mi "2Mp F ' ma " 2M P F ' m3 ~ 2Mp~ F ' 

5 _ 2g 2 F 5 1 

m ~ SM^M 4 ]^] m 

The last two of the above equations can be combined to eliminate (#23) and yield 

_ 5 vV F 5 v 

m 384tt 2 M 4 p^m z M P F' 1 ' 

a form that will be used in our numerical analysis. 

It is remarkable that the angle Qyi does not arise in Eq. (30). In the left hand sides of 
Eqs. (24) — (29) we have three independent neutrino mass eigenvalues and there are three 
independent parameters on the corresponding right hand sides. The three parameters can 
be expressed as: 

v - l(Xl2)l v - 1(3:23)1 v - ^ (32) 
Xi - — 7^- > X2 — — 7^-, X3 — ~ = 777- (o^J 



Upon using the relation (xi 2 ) = —(£13) in Eqs. (24) — (29), we are left with four equations. 
Any three of them can be used to solve \\-> X2 and X3 i n terms of m 1; m 2 , m 3 , c and s. On 
substituting the values of the x' s i n the fourth equation, we obtain a constraint equation 
among mi, m 2 , m 3 , c and s. The latter is automatically satisfied for mi = m 2 irrespective 
of the value of #i 2 . Thus the near-equality of two mass eigenvalues, basically reflecting 
the smallness of the mass splitting required by the solar neutrino deficit (as compared to 
that necessitated by the atmospheric neutrino shortfall), causes #i 2 to disappear from the 
solutions. 

The above feature can perhaps be motivated by symmetries of the neutrino mass matrix. 
As has been noted in recent works [T7j, when one sets 6* 23 = 7r/4, 6>i 3 = 0, and further neglects 
the mass splitting m 2 — mi, then the mass matrix becomes invariant under the successive 
interchange of the second and third rows, and the second and third columns. This symmetry 
is found to be independent of the value of # 12 : a feature to which the observations of the 
previous paragraph can be related. 
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Figure 2: /x — >■ ej, an FCNC process is shown at one loop level. 



3.4 Constraint from /i — > e"f 

Special care has been taken in our formulation to ensure a diagonal form for Fx tj so that 
FCNC processes are suppressed. However, a strongly constrained process like the radiative 
leptonic decay fx — > e^ can still receive one-loop contributions via two insertions of the L- 
violating sneutrino mass, as shown in Fig. 2. In order to estimate such a contribution, one 
can compute the corresponding amplitude: 

16tt 2 x 30m VM P ; ^ ; V 2 / v ^ w ' v ; 

where p, p' are momenta of the incoming muon and the outgoing electron in Fig. 2 and 
q = p' — p. Here V^ e \ are summed mixing matrix elements that enter the corresponding 
chargino-lepton-sneutrino vertices. In the above equation we have used A ~ Mew ~ m - 
By comparing this expression with, say the Standard Model amplitude for the /i — > e r y 
transition, we notice that there is an additional suppression factor of (J^j ~ 1CT 7 . This 
factor is small enough to automatically ensure a sufficiently low rate for the process /i — > ej. 



4 Different scenarios of neutrino mass hierarchy 

From the four equations (30), we notice that low-energy observables are ultimately controlled 
by three parameters of the model, namely, | (^12) I ? 1(^23)1 a- n d F. In order to fix them (or 
the ranges they lie in), though, one needs to know the values of the neutrino masses, along 
with the mixing angles. However, while the mixing angles are experimentally known to be in 
certain allowed ranges, all that we can claim to know so far about the masses are the mass- 
squared differences Am^ 2 and |Am| 3 |, corresponding to the solar and atmospheric neutrino 
deficits respectively. Their allowed ranges of values, together with those of the mixing angles, 
can be found, for example, in [ID] . Based on these ranges, all three scenarios, namely, normal 
hierarchy, inverted hierarchy and degenerate masses [IBJ, can be constructed in our model. 
Each of these places the individual mass eigenvalues within specified ranges. On using them, 
one can obtain the allowed ranges of the model parameters mentioned above. In the process, 
simultaneous use can be made of the fact that the SUSY breaking mass parameter fh is 
bounded from above if there are observable TeV-scale superparticles. Similarly, a lower 
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bound on rh can be imposed from negative superparticle searches with present accelerator 
data. Using the expression for m, the allowed space for the VEVs of the components of the 
gauge singlet chiral superfields can be further constrained. Thus one can check whether the 
VEVs of the scalar and auxiliary components of the X^, as required by neutrino masses, are 
consistent with the expected scale of SUSY breaking and the value of the /z-parameter. The 
self- consistency of the entire scheme gets established in this way. 

For illustration, we take the lower and upper bounds on rh to be about 100 GeV and 2 
TeV respectively. The allowed mass-squared difference ranges from oscillation data are 

Am 2 21 = (8.0 ± 0.3) x 10~ 5 eV 2 
|Am 2 2 | = (2.5 ±0.3) x 10" 3 eV 2 . (34) 

We have taken Mp = 2 x 10 18 GeV in our numerical analysis. The results presented below 
show the minimum value of \/F to be above 5 x 10 9 Gev, corresponding to the lower limit 
on m, which has its justification in the Large Electron Positron (LEP) collider results. 



4.1 Normal hierarchy 

This scenario corresponds to 

m! ps m 2 ~ \J Am 21 , m 3 ~ sJ\Krn^\. (35) 

In figures 3(a), (b) we have shown the allowed regions in the | (^12) | — VF and | (^23) | — ^J~F 
planes corresponding to the 3a range of ^Am^ as well as of \J\Am 32 \. On using the lower 
and upper limits of m, mentioned earlier, \/F is found to range between k 5x 10 9 GeV 
and 5 x 10 10 GeV. The scalar VEVs, on the other hand, are found to lie in the range of 10 11 
— 10 12 GeV. Finally, in figure 3(c) we have plotted rh against \fF using Eq. (31) for the 
lower and upper limits of m 3 at the 3a level. If F has to be related to the SUSY breaking 
mass terms, it is desirable to have it in the high side of the allowed region shown here. Thus 
values of F like a few times 10 10 GeV, and therefore rh somewhat on the higher side of the 
permissible range, are therefore favored in this model, given the accelerator search limits on 
superparticles. 

The next point to note is that while m 2 and m 3 have specific lower as well as upper limits 
in the normal hierarchy scenario, mi could, in principle, go down to zero. Nevertheless, the 
difference between mi and m 2 being quite small, the allowed region is restricted to be so 
narrow that it can be almost called fine-tuned. As discussed below, the situation is somewhat 
different in this respect in the case of inverted hierarchy. 

4.2 Inverted hierarchy 

Here we have 

mi ps m 2 ~ J\Amj 2 \, m 3 < J\Aml 2 \. (36) 
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Figure 3: Normal hierarchy: the top left (right) panel shows the allowed region in the 
| (£12) | (| (^23) |) — VF plane. The center lower panel shows the dependence of the SUSY mass 
scale fh on y~F in the allowed region of 777.3. 
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Figure 4: Inverted hierarchy: the left panel shows the allowed region in the |(#23)| ~ V F 
plane. The right panel shows the variation of the SUSY mass scale m against y/F in the 
allowed range of 777.3 ■ 
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(x 23 )| replaced by |(xi 2 )| in the y-axis. In addition 
Since we know that 
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Am§ 2 | 



gives us the same plot as figure 3(b), with 
£23)1 is plotted against \/~F in figure 
0.05 eV, we have allowed a maximum of 7723 = 0.01 eV. 



The minimum value of 777,3, on the other hand, could in principle be zero. In this case, too, 
for each m 3 value there is a limiting upper bound on F coming because 777 < 2 TeV. The 
corresponding parameter range is represented by the shaded area in figure 4(b), where 777 is 
plotted against \fF using Eq. (31) upto a maximum value of 7773 = 0.01 eV starting from 
7773 = 0. The interesting point to note here is that 7773 has no specified lower limit in this 
scenario. As a result, the allowed regions in the parameter space are much wider and less 
fine-tuned compared to the normal hierarchy scenario. Our conclusion, therefore, is that the 
inverted hierarchy scenario allows a larger flexibility of high-scale parameter combinations 
in the scheme adopted here. 



4.3 Degenerate masses 

This case corresponds to 

777! Hi 777 2 ~ 777 3 3> ^J[Am^\, (37) 

leading to | (a; 12) | ~ | (#23) | along with the requirement the actual masses shoul be significantly 
greater than the mass-square separations. Since \J\ Am| 2 | ~ 0.05 eV, we have to take 7773 ^> 
0.05 eV. On the other hand, since these are Majorana neutrinos, there is an upper bound 
of about 0.3 eV on the lightest mass from neutrinoless double beta decay as well as from 
cosmological constraints JUj- Therefore, we have plotted |(s 12 )| against yF for m 2 ranging 
from 0.1 eV to 777 2 = 0.3 eV, showing the allowed range as shaded area. This is shown in 
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Figure 5: Degenerate masses: the left panel shows the allowed parameter space in the 
I(^i2)| — VF plane. The right panel shows the variation of the SUSY mass scale rh against 
\F in the allowed range of 777,3. 



figure 5(a) with the usual constraints on F coming from rh. In figure 5(b), the allowed region 
in the rh — \[F plane is shown with m 3 ranging between 0.1 eV and 0.3 eV. 



4.4 Overall observations 

After analyzing these three cases, we can compare their impact on the parameters of our 
proposed model vis-a-vis the same on other similar models put forward in the literature. In 
the model considered in [3], for example, there is effectively one parameter which is F « (x) 2 . 
This is because just one right-chiral neutrino was considered there. That is why, despite the 
occurrence of both seesaw and radiative masses, the former are negligibly small in magnitude 
as compared to the latter. We have a more general (and natural) picture with three right- 
chiral neutrino superfields. There are consequently three unconnected parameters F, |(^i2)| 
and I (^23) I- These give us more freedom enabling us to treat the seesaw and the radiative 
masses on the same footing. The plots shown above for the three different neutrino mass 
cases imply that the different parameters which enter are in the expected range, thereby 
demonstrating the self-consistency of our model. 

It should be noted that the SUSY breaking scale rh is related to F by the relation 
rh ~ -j^, and that they are not entirely independent. The results presented in figures 3—5 
confirm that such a dependence is consistent with the requirement of neutrino masses and 
mixing. Some additional constraints may be required; for example, for rh on the higher side, 
one may be restricted to relatively larger values of F. However, it is impossible to be more 
exact in the absence of precise knowledge of the coupling strengths and other numerical 
factors in the hidden sector. Naively, eq. (31) is consistent with rh ~ for m.3 in the range 
.01 — .1 eV. This, in principle, restricts the inverted hierarchy scenario a little bit, although 
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it is difficult to be very precise, for reasons already mentioned. 

It is clear from the expression of Eq. (30) for m that our model connects very light 
neutrinos to TeV-scale massive particles. The latter include superparticles such as neutrali- 
nos and sneutrinos as well as right-chiral neutrinos. Further experimental information on 
neutrino masses, specifically the fixation of the hierarchy scenario, will therefore enable us 
to indirectly probe such yet undiscovered particles. At the same time, we have guidelines 
concerning the SUSY breaking sectors, especially the non-renormalizable terms that may 
have other ramifications such as explaining the worrisome /i-problem. 

5 Other possibilites 

In this section we briefly comment on two other representative scenarios where nonrenormal- 
izable interactions may be invoked to explain the observed pattern in the neutrino sector. 
We emphasize, however, that none of these has any bearing on the conclusions presented in 
the last two sections. We include these remarks mainly for the sake of completeness. 

In §3.2 we presented a superpotential which led us to diagonal A-terms and the consequent 
suppression of FCNC effects. This requires nonvanishing F-component VEVs only for the 
diagonal elements of the array of superfields Xy; the nondiagonal members of the array 
should have VEVs of the scalar components only, in order to generate off-diagonal elements 
of the neutrino mass matrix. One may be curious to ask whether, for the diagonal elements 
Xa, one can have nonvanishing VEVs for both the scalar and the auxiliary components and 
if such be the case, what their implications should be. 

There are some models based on Polonyi fields in which SUSY breaking has been achieved 
with a vanishing cosmological constant ^H]. In these supergravity-inspired models, on 
integrating out some additional chiral quark fields at a scale A, an effective superpotential 

W = XA 2 Z (38) 

is obtained at a scale below A. Here A is a constant 0(1) and Z is a Polonyi superfield. 
The scalar potential constructed therefrom yields (Fz) = AA 2 . At the same time, super- 
gravity effects lift the flatness of the direction (z) = so that both the scalar and auxiliary 
components of the chiral superfield Z have nonvanishing VEVs. 

Though one has in the past appealed to such models 13 EI, they do pose some 
difficulties in our case. A superpotential of the above kind cannot be used for off-diagonal 
elements of Xy, since those would generate large A^s for i ^ j and threaten to enhance 
FCNC rates. Thus the superpotentials for the diagonal and off-diagonal members would 
look very different in such a case, thereby raising doubts about the legitimacy of using the 
components of Xij as fields of a similar type. What we have done, on the other hand, does 
not raise such questions, since the complementarity of and Fx it is decided essentially 
by the relative magnitudes of two sets of mass parameters (/i, //) of the same order, where 
some fluctuation is quite natural. 

Another possibility [3] lies in considering the effective SUSY Lagrangian 

C eff = [XNN + LNH U } F + h.c, (39) 
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having one right-chiral neutrino to explain tiny neutrino masses. Eq. (39) can be readily 
generalized to include three right-chiral neutrinos, i.e. 

C eff = [XijrfW + XAjUWH^p + h.c, (40) 

where there is a summation over the flavor indices i,j and A« are Yukawa couplings. The 
Kronecker 5 in the second term is to suppress FCNC processes. This Lagrangian can be 
justified on the basis of the global symmetry Gp x G, where Gp = SO(3)p and G = 
U(1)r x U(1)l- For SUSY breaking, consider the hidden sector superpotential 

W = £ Wy, (41) 

where 

W i3 = .S; ;! .V, ; .V, ; - 4) + .YjV, ( + X%Y ir (42) 

The charges for various fields under G in this case are: 

^•(2,0), ^(0,2), ^(0,-2), 3^(2,4), (43) 
^■(2,-4), A^(l,-1), i?(l,l), if tt (0,0), H d (0,0). 

The minimization of the scalar potential, arising from this superpotential, yields (x^) ^ 
for all After SUSY breaking, Dirac and right-handed neutrino masses are generated as 

i m D}ij = ^ A *<%' l m Rh = ( X ij)- ( 44 ) 

Hence the seesaw mass is 

m v = —m D m~^m r D . (45) 

For (xij) ~ 10 11 GeV and Aj ~ 0.01, we obtain neutrino masses m u ~ 0.1 eV. Following the 
analysis given in the earlier sections, we can explain the observed bilarge neutrino mixing as 
well as different hierarchies in some parameter space of (x^). However, no connection can 
be made in this scenario between the tiny neutrino masses and the TeV-scale particles of 
MSSM. This is since only (x) enters the game and no (F). Moreover, a global symmetry of 
the form U(1)r x U(1)l does not allow the triple-X higher order term which can generate 
the AL = 2 sneutrino mass, cf. Eq. (9), in this particular model. Thus there can be no 
radiative contribution, at least not in the lowest orders. So one is unable to use here the 
full potential of such a scenario which, in our case, has meant a considerable widening of 
the parameter space through an interplay of seesaw and radiative effects, making our model 
more accommodating and natural. 



6 Summary and conclusion 

With a broken supersymmetric theory, we have considered a general scenario where the 
neutrino mass matrix is constructed through a combination of the seesaw mechanism and 
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radiative effects. All the agents behind the mechanism for this generation come from (a) 
a sector containing nonrenormalizable interactions in the superpotential and (b) lepton- 
number violating terms. We have used an array of gauge singlet chiral superfields for 
this purpose. In terms of these, we have constructed a superpotential which allows both the 
above types of contributions, while ensuring FCNC suppression. Right-chiral sneutrinos are 
found to have as much of a role in the process as the corresponding neutrinos. It should 
be noted that the angle 6 13 vanishes on using the bilarge mixing matrix in equation 2. A 
small but nonvanishing value of #13 requires one to modify equations (24) — (29), although no 
quatitative change in the conclusions is expected. We have then taken in turn the neutrino 
mass scenarios of normal hierarchy, inverted hierarchy and degenerate neutrinos. Using the 
masses answering to each scenario, we have traced out the allowed region of the parameter 
space of the involved high-scale physics, given in terms of the relevant VEVs of the scalar 
and auxiliary components of the superfields X^-. Numerically, these are seen to allow a self- 
consistent region of the parameter space. While the scenario of normal hierarchy (and, partly, 
that of degenerate neutrinos) forces us into somewhat fine-tuned zones of the parameter 
space, the inverted hierarchy picture allows a considerably larger region. With forthcoming 
laboratory measurements and cosmological observations hopefully deciding among the above 
mass patterns, connecting experimental observables to high-scale physics may become a 
realistic proposition, especially if Nature indeed proves to be supersymmetric. 
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